Normal- Superfluid Interface Scattering For Polarized Fermion Gases 
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We argue that, for the recent experiments with imbalanced fermion gases, a temperature difference may 
occur between the normal (N) and the gapped superfluid (SF) phase. Using the mean-field formalism, we 
study particle scattering off the N-SF interface from the deep BCS to the unitary regime. We show that 
the thermal conductivity across the interface drops exponentially fast with increasing h/kgT, where h is 
the chemical potential imbalance. This implies a blocking of thermal equilibration between the N and the 
SF phase. We also provide a possible mechanism for the creation of gap oscillations (FFLO-like states) as 
seen in recent studies on these systems. 

PACS numbers: 03.75. Ss, 03.75.Hh 
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Introduction — An electron approaching a normal- 
superconducting (N-SC) interface from the normal side 
with energy E < A, where A is the superconducting 
gap, has insufficient energy to excite quasiparticles inside 
the SC and is therefore reflected. In a SC, the relation 
A <^ Ep (where Ep is the Fermi energy) constrains the 
momentum transfer during interfacial scattering to be much 
lower than 2hkp, required for normal (specular) reflection. 
Consequently, Andreev reflection occurs: the electron pairs 
with another electron of opposite momentum, forming a 
Cooper pair. A "reflected" hole is left on the N side, which 
follows the time-reversed path of the incident electron. An- 
dreev used this process to describe the unusual heat con- 
ductivity found in N-SC junctions and in the intermediate 
state of superconductors [1]. To date, the transport prop- 
erties of such structures remain the subject of intense re- 
search 12, S]. On the other hand, recent experiments 10. S] 
have probed superfluidity in ultracold fermionic mixtures, 
where the possibility arises of having different chemical 
potentials for each species and controllable interspecies in- 
teractions. An accurate theoretical prediction of the ob- 
served density profiles is still not available, despite intense 
theoretical activity ||E 0, S i, [H [H HI. Whereas the 
incorporation of a normal-superfluid (N-SF) interface ten- 
sion may settle this question [4, Jj\ for the Rice experi- 
ments, the interpretation of the MIT experiments remains 
unclear. A recently proposed theory (E] appears to explain 
all observed features, although the temperatures required 
for agreement with data from the MIT experiments are five 
times higher than observed. 

In the following, we argue that the presence of a N-SF in- 
terface is likely to block thermal equilibration for these ex- 
periments, possibly inducing a temperature difference be- 
tween the two phases. Incorporating such a temperature 
difference in the existing models may provide the key to a 
complete understanding of the experiments. 

In this Letter we focus on two question, to wit: 1) What 
are the possible reflection and transmission mechanisms at 
a N-SF interface? The main features distinguishing this 
system from a SC are the difference in the chemical po- 



tentials and the variation of the gap from A ^ Ep to 
A ~ Ep as one tunes the interactions from deep BCS 
to the unitary regime. 2) What is the relevance of these 
findings for recent experimental and theoretical works? 

We find a rich variety of interfacial scattering processes, 
depending on the energy and perpendicular momentum of 
the incoming particle. The scattering of the particles off 
the interface results in a striking decrease in the thermal 
conductivity for temperatures T < 0.05 Tp (at unitarity). 
We also argue that upon Andreev reflection, particles and 
holes interfere so as to cause gap oscillations near the inter- 
face, which are observed in recent numerical studies [12]. 
We postpone the description of calculational details and the 
case of unequal fermion masses jlSh . 

The system under study consists of two fermionic 
species a and b with equal masses m and chemical poten- 
tials for i = a, b, trapped by a potential V(r) ; this gives 
rise to an effective chemical potential /Ui(r) = — V(r) 
(henceforth denoted by /Xj). The Bogoliubov-de Gennes 
(BdG) or Blonder-Tinkham-Klapwijk equations ^ [13] 
give a satisfactory description of this system. Denoting 
the particle-like and hole-like eigenfunctions of species 
i = a, bby k and Vi,^ respectively and following the 
standard procedure, one finds: 
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withTia.b = —V^/2m — /i=F^ where /i = (fia+fJ-b)/'^ > 
0, h = {fia — fJ'b)/'^ > and we set h = ks = 1- One 
obtains the second set of BdG Eqs. by interchanging a and 
b. As seen from Eqn. ([Tjl, incoming a-particles are cou- 
pled to 6-holes and vice versa. Now, the prerequisite for 
an interface to be present between a SF with symmetrical 
a and b densities and an asymmetrical N phase is the ex- 
istence of a first-order transition between the two phases. 
In the experiments, if a SF is found at the trap center, /j, 
decreases upon approaching the trap boundary and probes 
the {fia, fJ-b) phase diagram, possibly inducing a crossing 
of the N-SF phase boundary. This happens when locally 
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FIG. 1: The N-SF interface (thick vertical line) with the a, b 
spectra on the N side and the gapped a, (3 spectra in the SF. The 
long dashed line cuts the spectra at particle-like (filled dots) and 
hole-like (empty dots) quasiparticle states, all having the same 
energy. An incoming a-particle (curly arrow) with momentum 
kp and energy _E > A — ft has up to four scattering channels: the 
Andreev reflected fc/j hole, the specularly reflected —kp particle 
and the transmitted hole-like — fc'' and particle-like fc^ states. 



there occurs a balance between the energy gained by creat- 
ing, on the one hand, a gap and on the other hand, a density 
difference [17]. At T = 0, this local coexistence condition 
is well-approximated by the Clogston limit [J 3,1 : 

A = \/2/i. (2) 

For fixed interaction parameter kpa, where a is the scatter- 
ing length, the first-order transition persists at finite temper- 
ature up to a tricritical point [^JJJ. The Rice experiments 
are indeed accurately described using a theory incorporat- 
ing a first-order transition S 0, 0, 'Sj ; the MIT results are 
also suggestive of a first-order transition. 

To model the N-SF interface, it is natural to start with a 
geometry wherein the x = plane separates the N from 
the SF region: 

A(r) = G(x)A, (3) 

with Q the Heaviside function. This choice allows one to 
find the general solutions of the BdG Eqs. ([B in the N and 
SF phase. They read: 
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Here, the sub- and superscripts n, s, p and h denote nor- 
mal, superfiuid, particle-like and hole -like respectively. We 
split the vectors k into their components parallel to the wall 
k|| = (0, ky, kz) and the x-component kp which relates to 
the A;^ '' by the BdG Eqs.: 
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FIG. 2: The various scattering regions for an a-particle incident 
on the interface from the N side, as a function of its energy E and 
= kp/2m (kp is its momentum along the i-axis). The heavily- 
shaded region is energetically forbidden, while complete reflec- 
tion occurs in the lightly-shaded regions. Particles in region VI 
may scatter to all states indicated in Fig.[T]by curly arrows. Above 
line 4 (regions IV and V), Andreev reflections do not occur, and 
above curve 3 (region IV), hole-like excitations are also impos- 
sible. The numbered curves are 1: = 2: ^p — fia + E, 
3:Cp = ?+and4:ep = 2(£; + ft). 



where (,±= E + h± ^J{E + hy - A^. The amplitudes 
U and V from solutions (|4]l are determined by matching the 
wavefunctions and their derivatives at x = [15]. Figure [T] 
shows the quasiparticle spectra in both the N region (left 
panel) and the SF region (right panel). In the former one 
recognizes two approximately linear energy branches a and 
h near the Fermi surfaces, while in the latter the a and fi 
spectra are gapped by A — /i and A + /i. Since we take 
> fJ-b, the states belonging to the a (lower) spectrum 
are the ones composed of a-particles and 5-holes. 

Scattering regimes — Consider an incoming a-particle 
from the N side with energy E and momentum k; it is indi- 
cated in Fig.[T]by a curly arrow and its motion is completely 
characterized by E and = k^/2m. The crossover of ei- 
ther of the fc^ '* or k|| from real to imaginary (or vice versa) 
signifies a change in the scattering mechanism. The regions 
in the E — ^p plane corresponding to the various such scat- 
tering processes are shown in Fig. |2l we now briefly de- 
scribe them. If £" < A — /i, the incoming particle has 
insufficient energy to excite quasiparticles inside the SF; 
it is completely reflected as a superposition of a particle 
and a hole. This region is labelled I in Fig. |2l The sit- 
uation where E < — /i^ is physically forbidden for 
obvious reasons; it is labelled II. Next, consider energies 
above the threshold A — h. For < ^_ (region III), 
the incoming particle can excite neither particle -like nor 
hole-like quasiparticles in the SF and must therefore un- 
dergo complete reflection, even though it appears to have 
sufficient energy for transmission. The situation is akin 
to that of a quantal particle incident on a potential step; 
if the angle of incidence 9 exceeds a critical angle, the 
particle is reflected. Here, the critical angle is 9^ with 



3 



tan Of = ^J{E + h + - 1, in which 6' < cor- 

responds to < ^- . For angles 9 satisfying 9^ < 9 < 6^ 
(region IV), or < < particle-like states may 
be excited in the SF, in contrast to hole-like states. When 
< 2{E + h) (regions IV and V), there are no reflected 
holes; that is, Andreev reflection does not occur. Thus, 
in region V, both particle-like and hole-like excitations are 
present in the SF, but Andreev reflection is impossible. Fi- 
nally, in region VI, both particle-like and hole-like excita- 
tions, as well as Andreev and normal reflection are allowed. 
To summarize, the lightly-shaded regions in Fig.|2]describe 
a-particles which undergo complete reflection, while the 
rightmost, heavily-shaded region is unphysical. Only in- 
coming particles in the unshaded regions may excite quasi- 
particles inside the SF. As for incoming holes, they are only 
transmitted in a region analogous to II in the diagram. Ar- 
guments very similar to the preceding, applied to incoming 
6-particles and a-holes (therefore the /3-channel), lead to a 
diagram identical to that of Fig.|2]but with h — > —h. 

The £,p-E diagrams of the deep BCS and the unitary 
regime have a different topology. In the deep BCS regime, 
the relation 2 A <^ A + /i implies that region VI of Fig. |2]is 
by far the most important. Quasiparticle reflections occur 
mostly via the Andreev mechanism which involves parti- 
cle, but not energy, transport across the interface. In con- 
trast, at unitarity where 2 A > A + (since A ~ 1.16/i), 
regime VI does not even exist. This means that quasiparti- 
cles with energy above the threshold A — h cannot undergo 
Andreev reflection but only normal (specular) reflection, in 
which neither particles nor energy are carried across the 
interface. 

Thermal conductivity k — In order to study the trans- 
port properties, we seek to relate the amplitudes U and V 
of Eqs. ^ to the transport coefficients, based on a conser- 
vation law for the current. From the BdG Eqs. it readily 
follows that the density ^^(r) = |tiaP + and the 
quasiparticle current 

2m 

satisfy the continuity equation. The transmission coeffi- 
cient S{E, ^p) of an incoming particle of energy E and 
momentum k is defined as the ratio of the transmitted to 
the incoming current along the x-axis. S vanishes in re- 
gions I, II and III whereas for energies slightly above the 
transmission threshold E ^ A — h (for regions V and VI), 

S{E, Q oc y'E-iA-h), (6) 

similarly to the case of particles scattering from a Hartree 
potential of height A — h. Next, one can write the heat flux 
through the interface and to the i = a, P branch as: 

^-^^ = £^Y1 I <iEEf{E)S{E,C,), (7) 

s—pji 

where f{E) is the Fermi-Dirac distribution and the inte- 
gration is performed over the ^p-E and S,h-E planes (with 
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FIG. 3: The thermal conductivity across the N-SF interface k di- 
vided by the normal-phase conductivity kjv against T/A for the 
unitary, BCS and deep BCS cases. For T < O.IA, k/kjv drops 
dramatically (notice the logarithmic scale). The curve represents 
the analytical result obtained by the use of the Andreev approx- 
imation, Eqn. At unitarity, A = OWTp, in the BCS case 
A = 0.25rir and in the deep BCS case A = 0.002 Tp. 



= kf/2m). In equilibrium, the N-SF flux is balanced 
by an equal SF-N flux. A small temperature bias on the 
N side will induce a net heat flow Q. By the Kapitza ap- 
proach, the latter is Q = k 5T with the heat conductivity 

AC = diW.,^^ + W^^p)/dT. 

We have analytically calculated the transmission coeffi- 
cients for all regions of Fig. |2l To fix the values of h, 
and Tp = kp/2m = (37r^n)^/'^/2m for given kpa (n 
is the density in the SF), we have used the gap and num- 
ber equations, as well as the exact T = coexistence 
condition |fl3ll . The resulting ratio of k to the conductiv- 
ity in the N phase, k^, is shown in Fig. [3] as a function 
of T/A. We find that k/k^ decreases drastically below 
T ~ O.IA 1 18]. Remarkably, h/k.^ is almost independent 
of the interaction parameter kpa. This can be understood 
as follows. For low enough temperature, T <^ A, only in- 
cident particles (and holes) with energies slightly above the 
threshold E A — h contribute to the heat conductivity 
because of their low statistical weight / oc e^^^^'^^^'^ ap- 
pearing in eqn. Furthermore, relation (O shows that S 
displays the same behavior in regions V and VI. The main 
effect of varying 1/kpa from large and negative (deep 
BCS) to zero (unitarity) is to move the boundary between 
regions V and VI to the right. Both the strong exponential 
decay with A /T and the square root dependence of S on 
E — (A — h) are unaffected by the variation of 1/kpa, 
hence the very similar behavior of k/kn for all regimes 
under study. 

These considerations are exemplified by the following 
calculation in the deep BCS regime, where the Andreev 
approximation may be used to obtain an analytical expres- 
sion for K. The approximation is to equate all the k^''/^ wave 
vectors and take the temperature to satisfy T ^ A ^ /i. 
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We obtain 
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which amounts to Andreev's result when /i = H]]. The 
summation denotes the excitations of the a (-) and /3 (+) 
states in the SF. The energy carried by the (3 branch is a 
factor e~^''/^ lower than that of the a branch and due to 
coexistence h A ^ T, it can be neglected. Equa- 
tion dSJ demonstrates the exponential decay of k with de- 
creasing temperature. Although the Andreev approxima- 
tion is invalid for the aforementioned experiments |0, IH], 
formula ([8]) as plotted in Fig. [3] provides a good estimate 
of K, even beyond the BCS regime lfl9h . 

FFLO states — Andreev reflection may also produce 
effects similar to Fulde-Ferrell-Larkin-Ovchinnikov states 
(FFLO) as seen in Ref.Hl'. For the case of a SC, McMillan 
calculated 1 16] the first corrections to the gap function (O 
due to self-consistency and found, at T = 0, a decaying os- 
cillation of A(x) on the N side. This stems from the phase 
difference ib(/Cp — k^) between the wavefunctions of in- 
coming particles (holes) with momentum kp [k^) and en- 
ergy E < A, and their Andreev reflected holes (particles) 
of momentum A;^ {kp)\ the wavefunctions of these deter- 
mine the gap profile by A{x) oc ^j^. Ua^\^{x)vl_]^{x). We 
speculate that this also happens beyond the BCS regime. 
In particular, the numerically observed T = FFLO states 
for trapped gases at unitarity may be a consequence of the 
presence of the N-SF interface |12]. Note that the FFLO 
state is not thermodynamically stable for a homogeneous 
system at unitarity [9]. 

Discussion — In a SC, the heat conductivity also has a 
lattice component, which dominates the electronic compo- 
nent at low temperatures; the absence of such a component 
in the system under study makes the decrease of the con- 
ductivity more significant. Addition of the Hartree-Fock 
potentials amounts to a change of the coexistence condi- 
tion dJll and a mere horizontal shift of lines 1 through 4 
in Fig. |2] thereby preserving the energy gap (region I) and 
thus the main conclusions of our paper. We expect that im- 
provements on our approximations, which are the use of 
the one-channel model, the gap profile (|3]l, and the temper- 
ature dependence which is solely contained in the Fermi- 
Dirac function, will affect the results only quantitatively 
and mostly at unitarity. Finally, note that the values we 
used for /i, /x and Ep are those at the interface. Estimat- 
ing that Tp at unitarity (see caption of Fig. [3]l is equal the 
reported iBll, one concludes that the drop of conduc- 
tivity sets in at T = 0.05 Tp, from which follows that our 
arguments indeed apply to the current experiments. 

Conclusion — Scattering of quasiparticles off the N-SF 
interface, as summarized in Fig.|2] gives rise to the follow- 
ing effects. The thermal conductivity across the interface 
decreases rapidly with decreasing temperature at the exper- 
imentally realized temperatures, and is described by equa- 
tion dH. This implies that thermalization slows down, and 



a temperature difference across the interface can appear; 
the incorporation of this could lead recent models to better 
agreement with current experiments. We call for such in- 
corporation, as well as for a separate temperature measure- 
ment of the N and SF phase since, as we have argued, they 
may not be the same. In addition, we argue that reflections 
of quasiparticles off the interface can induce modifications 
to (BCS regime), or even cause (at unitarity) gap oscilla- 
tions, usually identified as FFLO state. 
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